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Abstract 



Using the statistical approach additional light partons have been introduced in the Intrinsic 
Charm model. Explicit expressions are obtained for charmed and light quark distributions. A 
comparison is made with standard IC model predictions. It is shown that such a modification 
of the IC model leads to better agreement with the pN and ttN data on J/tp production at low 
energies and allows to avoid a contradiction to experimental data at high energy. 
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1 Preface 



Evidence for an excess of charm production at high xp observed in a number of experiments 
on proton-proton |l| and muon-proton ||] collisions stimulated the development of the model 
with an assumption that the hadron wave function decomposition may contain the | uudcc > Fock 
state component. The first theoretical estimation for the charmed particles production due to 
nonpertubative charm component inside the proton has been made by S.Brodsky et al. in |3] using 
old-fashioned pertubation theory. The standard Intrinsic Charm (IC) model predicts very hard 
distributions for charmed particles due to the fact that the charmed quarks in the proton carries 
a large fraction of the proton longitudinal momenta. But some last experimental data [[h]] on the 
J/ip production in pN interactions at P$ = 800 GeV/c show that the distribution predicted by 
the model § strictly contradicts to the experiment at large xf- 

If the above assumption about excitation of the higher Fock states with heavy quarks is valid there 
must be much larger probability to create additional light partons. Therefore we have to consider 
also the Fock states with arbitrary number of light quarks and gluons. With natural restriction 
that this parton system as whole has to carry the hadron quantum numbers and sum of the parton 
momenta is to be equal to the hadron momentum. In this case the additional light partons restrict 
the available phase space and lead to effective decrease of the average longitudinal momentum 
of heavy quarks. Despite the recreated interest to IC problem no attempts have been made so 
far to consider this effect. This step has obvious physical basis and, as we show, leads to visible 
experimental consequences due to softer momentum distributions of the partons. The increase of 
the number of light quarks leads also to larger absolute value of total electric charge of the partonic 
system and to the larger cross section in the electroproduction processes. 

In present paper we exploit the statistical approach Q , derive expressions for distributions of heavy 
and light partons within framework of the non-covariant pertubation theory and compare the xp 
distributions of the J/ip particle to experimental data in irN and pN collisions. 

We shall speak mostly about the proton. For meson all the calculations are completely similar. The 
only difference is that in the meson case the lowest Fock state contains two light quarks instead of 
three ones and, therefore, number of integrations in the following expressions will be lesser by unit 
in comparison to the proton case. 

2 Description of the model 
2.1 The original IC model 

Further in present paper we shall denote as a Fock state \uudcc + n > containing one (cc) 
pair, three (two for meson) "valence" quarks, earring the hadron quantum numbers and n light 
partons, gluons and quark-antiquark pairs, with zero summary quantum numbers. In papers || it 
was proposed that the probability to create a Fock-state in the proton is being described by 
the expression: 

dW^( Xl ,x 2 ,...,x n+5 ) ~(m 2 -Y. — ) Ml - E ^ II dx ^ C 1 ) 

V i=i x% ) \ i=\ ) i=i 
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where the m± j are the effective transverse masses of the partons m ^ = mj + < p^ > , the < p\ > 
is the average transverse momentum squared, M is the proton mass and n denotes the number of 
additional light partons. 

If to regard valence quarks only, i.e., the Fock state = \uudcc > and to neglect mass of 
the proton M and all the mj_ for light quarks one easily obtains the one-particle distributions, 
normalized on unit: 



cf\x) = 1800x 2 ±(l-x)(l + 10x + x 2 ) + 2x(l-x) In x 



= 6(1-X) 5 

for charmed and light quarks, respectively. 



(2) 



Because we are interested only in the hadron Fock states containing two charmed quarks we can put 
in the folowing calculations the overall normalization Njq equal to unit for the sake of simplicity. 
To recalculate final results it is enough to multiply them by factor Nj^, extracted either from an 
experiment or from a theoretical model. This normalization is order of Njq ~ (0.5 — 1)% (see 

§§§!)• 



2.2 Arbitrary number of light partons 



In this section we use the statistical approach Q to introduce arbitrary number of additional light 
partons. Remind that we consider all the partons are being on the mass shell and use the non- 
covariant pertubation theory ||]. It is also necessary to keep in mind that the results obtained 
in present paper are, generally speaking, valid only for partons with sufficiently large momentum 
fractions ("valence" partons). The wee partons are being created due to bremsstrahlung and 
have different distribution. In present paper we don't consider their contribution and we use the 
condition J2i I dx fi(x) = 1. If we introduce additional sea partons following Eq.(|l]) and negelect 
transverse masses of light partons all the light quarks and gluons will have the same momentum 
distribution. Thus we can not distinguish the "valence" and the "sea" light quarks. The validity 
of the approximation (|l]) for many partons we consider later in this paper. 

The probability to create a final state is proportional to d^p/(2Trh) 3 — the number of cells in the 
momentum space. Because we don't consider transverse movement we can use the longitudinal 
phase space only. We define pk = Pmax x k, where P m ax is the maximal longitudinal momentum of 
the parton. Thus we can write for the probability to create either one gluon (g) or one quark (q) 
(remind, that fi = 1): 

dW ~ -^-Wg t qP max dx g = ag tq dxg tq , (3) 

Z 7T 

where, Wk is the statistical weight of the parton depending on its colour and spin. To simplify 
further formulae we introduced the entity a gA = (b/2ir) w g ^ q P ma x- 

For the on-mass-shell gluons we have eight colours and two polarization states. Similarly for each 
quark we have three colours and two polarization states. Thus we have: 
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Wg = 16, W q = 6. 



(4) 



If we neglect any dynamical effects we can put b g ~ b q = b ~ const. The parameter 6 is unknown 
and we at present have no indication on its exact value. But we can expect that it is order of unit. 
This parameter has to be adjusted either from experiment or from a theoretical model. E.g., in 
the bag model we are able, in principle, to calculate the necessary probabilities and distributions 
but it is completely unclear how will these values be transformed under the Lorentz boost to the 
infinite momentum system. Therefore the most direct way is a comparison to experimental data. 

We don't distingush the flavours of the (qq) pairs, i.e., we consider the S'C/(3)-symmetric sea. This 
suggestion does not play any role in our case. So we can write the expression for the probability to 
create a state with one (cc) pair, three quarks uud (two quarks for a meson), n g gluons and 
n q quark-antiquark pairs: 

n n a n 2n q Ji „2 n+5 n+5 

Here n = n g + 2n q , indecies "1" and "2" denote the c and c quarks, respectively. The factors 
l/n g \ and \/{2n q )\ take into account indistingushibility of the gluons and the quarks. 

To obtain the total probability for all the Fock states containing two charmed quarks one needs to 
perform summation over n, < n < oo, and the integration over all X{\ 

oo . 

1 = Z-\a g ,a q )Y, / dWi n \ Xl ,...,x n+5 ), (6) 

n=0 J 

where, Z~ l (a g ,a q ) is the normalization factor (we put the overall normalization equal to unit). If 
we perform the integration over all X{ in Eq.(||) we obtain: 

v ' n g +2n q =n f v 1' 

For the normalization factor we have: 



oo -, _rig 2n q 



n=0 n=0 \ n ^°>- n g +2n q =n n 9- ^l)- 

The summation over n = n g + 2n q can be carried out using properties of the binomial sums. As 
result we obtain for the proton: 

I 00 a n _|_ a n_ l 

Z p {a g ,a q ) = — t 6 T, = -5* 6 (a+,a_); a± = a g ± a,. (9) 

n=0 » ' 



Similarly for meson we have 
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Zm(o 9 , a q ) 



-E 

in i 



al + a r _ 



10 n!(n + 5)! 

n=0 v ' 



- S 5 (a + ,a. 
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where we denoted 



al + a n _ 



2 ^— ' n!(n + m) 

n=0 v ' 



E 



- |a + m/2 7 m (2^/a7) + o_ m/2 / m (2^/aT)| 



(10) 



Here, I m (z) is the modified Bessel function of the m-th order. 

In the considered model in difference to usual parton picture the distributions of the additional 
light partons does not look like bremsstrahlung spectrum and its value is finite at x — > 0. It follows 
from Eq.(P) that the coefficients in the sum have a maximum at some number n = n e ff. Thus 
we can speak about effective or average number of light partons in the hadron. Fig.|| presents the 
average number and average momentum fraction of light partons versus the P ma x m the proton at 
the above choice of the relative weights w g , w q . 
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Figure 1: The average momentum fraction (a) and the average number (b) of light partons versus 
maximal momentum P ma x for the proton. 

For light parton we obtain the following 1-particle distribution (for the proton case): 



qi(x;a) = - Z p 1 (a g , a q ) (1 - xf S 5 [a+(l - x), a_(l - x) 



(11) 



And for charmed quarks we have the 2-particle distribution: 



c 2 (x c ,x s ;a g ,a q ) = Z l {a g ,a q ) 



(l-x c -x 5 ) 2 S 2 [a+(l - x c - Xc),a-(1 - x c - x 5 )] , (12) 



One-particle distribution for charmed quark looks as 



ci(x c ;a g ,a q ) 



Z-\a g , a q ) x 2 c {l-x c f f 1 ^ V f ( \ l ^ S 2 [o+(l - x c )(l - y), o_(l - x c ){\ - y)] 



o [x c + (1 - cc c )yp 



(13) 



b s = 1.0 W 9 = 16.0 
b = 1.0 W= 6.0 



1 - No sea 

2- P m „= 10.0 GeV 

3- P ma ,= 100.0 GeV 




b s = '.0 W s = 16.0 
b= '.0 W= 6.0 



1 - No sea 

2- P mo ,= 10.0 GeV 

3- P mo ,= 100.0 GeV 
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Figure 2: Distributions for light (a) and charmed (b) quarks over x in the proton from Eqs.(ll), 
(13) at 

Pmax — 10) and 100 GeV. The solid line is the distributions predicted by the IC model 
without additional sea partons. 



Additional light partons carry zero quantum numbers and its number n is not invariant, i.e., 
depends on the reference frame. From the definition (0) and Figjl] we see that the average number 
of the light partons depends on the reaction energy. The increase in the energy of the reaction 
leads to the increase in the effective number of sea partons and to softer longitudinal distribution. 



We have to note that the growth of the < n > is faster than usual logarithmic law due to exponential 
growth of the number of cells in the momentum space. We assumed the constant distribution over 
the momentum for the unconstrained probabilities to produce a light parton. In general, it is not 
so and the parameter "b" can depend on the parton momentum. 



2.3 Self— consistency of the model 



Now we have to make a rough check of the self-consistency of the Eq.(|l|). Let us consider one 
Fock state with fixed number n of additional light partons. Perfoming integration in Eq.((l]) over 
necessary parton momenta it is easy to obtain the 2-particle distribution for cc pair: 

4 n) (*c,xe) ~ 7 X l X % 2 (l-x c -x 5 f+ 2 , (14) 

\X C "T x c j 

and respectively (n + 3) particle distribution for light partons: 
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g„+ 3 (xi,a?2 ! • • -,X n +3) 



These two expressions lead to 1-particle distributions: 



' n+3 



i=l 



c\ '{x c ) ~ x c (1 - x c ) / dy ■ 



o [x c + (1 - x c )yX 



for charmed quarks and 



for light quarks. 



q{ n \x) ~ (1-s) 



71+5 



(15) 



(16) 



(17) 



Average fractions of the proton momentum carried by a quark is: 

< x c > = ^nf, for c quark 

< x q > = for light parton. 



(18) 



All the above formulae have been derived in suggestion that effective transverse masses of heavy 
quarks are much larger than for light ones. It is really so but for validity of the approach we need 
to satisfy more restrictive condition for each Fock state F^: 



2 < m\ c > 
< x c > 



M 



2 y < m±i > 

fri < Xi > 



We take: 



2mi 



4.5 GeV 2 ; < x c > w 2 < x q > 



For < m\ q > < p\ q > 0.1 GeV 2 and for < x 9 > l/(n+7) one can obtain that for n < 12 the 
above approximation is valid within 30% accuracy. We hope that general formulae from previous 
subsection presenting the sum over all states work due to double factorial in the denominators of 
the expansion coefficients, which suppresses the contributions of the higher states. Of course, at 
small x q where the distribution q\(x) 3> c\{x) the validity of the approximation will be defmetely 
brocken but this model considers large values of x. 



So we conclude that our expressions ax6 self^consistent for sufficiently large Pmax 

< 20 (see 

Fig||(b)). 
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3 Application to J/ip Production 

3.1 The recombination model 



To estimate the longitudinal distributions of the J/ tp particles we use the recombination model 
In this model the differential cross section of final particle can be written as: 



da 

dXF 



/ dx\dx 2 F(xi,x 2 ) R(xi,X2]Xf), (19) 
Jo 



where, F{x\,x 2 ) is the two-particle distribution for quarks "1" and "2" with momentum fractions 
x\ and x 2 , respectively; R(xi,x 2 ;xp) is the recombination function describing the probability for 
two quarks to coalesce in the final meson with momentum fraction xp. In the simplest case: 

R(xi,x 2 ;xf) = S(xf — x\ — x 2 ), 

ensuring the longitudinal momentum conservation. In principle for charmed particles the primordial 
transverse momenta of the initial c quarks can reach large values (> 1 GeV) and have to be taken 
into account by obvious way (see, e.g., ||). But for the aim of present paper it is not important. As 
a result from Eqs.(|i~9|) and (|l~2^) one obtains for the xp distribution of J/i/j particles the following 
expressions: 

MJ i^ (pp) = ±Z- 1 (a g ,a g )x F (l-XF) 2 S 2 [a + (l-XF),a-(l-XF)], 

(20) 

^iS^ = m Z lh a 9i a q) x F( l - x F )S 2 [a + (l- x F ),a-(l- x F )}, 



for pp and irp collisions, respectively. For hadronic interactions we have in the center-of-mass 
system P max = yfs/2. If we normalize our cross section on specific value we have the only free 
parameter — the constant b from Eq.(||) At the moment we have no other choice as to put it equal 
to unit, 6 = 1. Note here that at larger values of b we have softer distributions for c quarks and 
for J/tp particle. 

In paper [|l(]] the J/ip production in the proton-nucleus interactions has been measured at -y/s = 
38.7 GeV and 0.30 < xp < 0.95. It is a region where we expect the contribution from the intrinsic 
charm. Based on estimations from [|llj and private communication to R.Vogt the authors of the 
paper [|l(]] made a conclusion that the IC model conradicts to their experimental data. They used 
a value for total cross section for IC J/ip production 0"t o t(J/V') ~ 1-8 nb. The present paper is not 
right place to discuss the absolute value of the J/ip total cross section. It is very complicated subject 
and needs special consideration. So we use the value 0.9 nb for a{J/ijj) in the forward hemisphere 
at FNAL energy. For lower energy data [Q| we normilized our expressions for measured total cross 
section in the forward direction for the so called " diffractive" component (see [l^]) of the J/ip yield. 

Results of this comparison are shown in FigsJH and § for b = 1. In these figures we also plotted the 
predictions of the standard IC model based on the same recombination model and with the same 
normalization. As one observes the predictions of the usual IC model can really be rejected for 
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above value of the total cross section at FNAL energy. On another hand predictions followed from 
our calculations don't contradict to the experimental data. At low energies both models agrees 
with the experimental data. 

This fact allows us to suggest that the contribution from intrinsic charm can be observed at suf- 
ficiently low initial energies and is practically negligible at high enrgy even at large xp. This 
comparison also confirms that parameter b ~ 1. 
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Figure 3: J/ip distributions over xp in (a) pp at Pq = 200 GeV/c and (b) in pCu collisions at 
P = 800 GeV/c. 
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Figure 4: J/ip distributions over xf in it p collisions at (a) Po = 150 GeV/c; (b) Pq = 200 GeV/c; 
(c) P = 280 GeV/c and (d) in vr+p collisions at P = 200 GeV/c. 
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